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<—< \ We obtain a semiclassical expression for the projector onto eigenfunctions by means of the Fred- 

holm theory. We express the projector in the coherent state basis, thus obtaining the semiclassical 
Husimi representation of the stadium eigenfunctions, which is written in terms of classical invariants: 
periodic points, their monodromy matrices and Maslov indices. 
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I. INTRODUCTION 

The precise test of semiclassical approximations in the presence of chaos is of grest interest to establish the limits 
of applicability of periodic orbit theory and its resummations. This test can be done in model systems both on 
approximations to the spectrum or to the stationary states. For the calculation of the spectrum the most efficient 
tool in this respect seems to be the spectral determinant jl]J^] and several calculations || have demonstrated that, 
given enough periodic orbits, the spectrum can be accurately represented semiclassically. However, a more sensitive 
test - and still a great challenge - is the semiclassical representation of single eigenfunctions. This includes the study 
of the scar phenomena [PI-ITo! and the eventual deviations from uniformity of eigenfunctions in accordance to the 
| Berry- Voros hypothesis [ |liyi2f| and Schnirelman's theorem [fj"3f . 

Just as for spectral problems, the use of Fredholm methods allows for the most efficient encoding of classical 
^f} " information in the calculation for single eigenfunctions |^,|l4|| . In this paper we review these methods and apply them 
f^) , to the calculation of Husimi distributions of stadium eigenfunctions. 

This paper is organized as follows. In Sec. || we review the Fredholm method for billiard eigenfunctions. Fredholm 
theory allows us to find the solution to certain type of integral or operator equations [pl| . For billiards these methods 



can be applied to the boundary integral equation. In Sec. Ill we make the semiclassical approximation that is 
based on the approximation of the traces and powers of the propagator as sums over the periodic points of the 
underlying classical system. The propagator itself is taken as Bogomolny's T operator Q. We choose the coherent 
state representation and obtain an expression for the semiclassical Husimi representation of the eigenfunctions in 
terms of classical invariants: periodic points, their monodromy matrices and Maslov indeces. In Sec. ^ we apply 
' this scheme for the stadium billiard. Our conclusions and perspectives are presented in Sec. |v|. 



II. FREDHOLM FORMULAE FOR EIGENFUNCTIONS 



Fredholm theory gives the solution to a certain class of integral equations, which can also be written as operator 
equations . A Fredholm integral equation of second type is 

x(q)=Xo(q) + X [ dq'T( q ',q)x(q'). (1) 



All the functions are defined in a finite domain. If the known functions Xo(o) an d T(q',q) are well behaved, the 
Fredholm alternative holds: there is a unique solution x with the same analytic properties or the homogeneous 
equation (xo = 0) has a solution. There is a set of complex parameters Xi for which the solution is not unique. In 
operator notation, the inverse of (1 — AT) exists if A ^ Xi. In this case, this inverse can be written as 

1 = M(A) 

1 - AT D(X) ' 1 ' 
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where the operator M(A) and the function D(X) are series in A. If T is a compact operator, D(X) and M(A) are 
entire in A and, thus, absolutely convergent. The explicit form for the series expansion in terms of powers of T is 
given below. In what follows we apply this general theory assuming T to be unitary and of finite dimension N. Both 
assumptions are justified in the semiclassical limit for the quantization of billiards [fij. 



A. Secular equation 

The fc's eigenvalues are given by the secular equation P{k) — det(l — T(fc)) = 0. We can expand this determinant 



as 

N 



P{k) = J2Mk), (3) 

71 = 

where the coefficients (3 n (k) are related to the traces of T(fc), b n (k) = trT™(fc), through 

1 ™ 

Pn(k) = --£&-,•(*)&,■(*)■ (4) 



n 
i=i 

Thus, knowledge of the traces up to a certain n max implies the knowledge of the coefficients (3 n up to the same n max . 
If T(fc) is unitary, P{k) is self reversive, meaning that its coefficients satisfy 

P N -j{k) = {—1) N Pj(k) detT(fc). (5) 

This condition alone forces the eigenvalues of T at fixed k to be symmetric with respect to the unit circle: if A is 
an eigenvalue, then 1/A is an eigenvalue too. Of course, if T is unitary, then this condition is automatically satisfied 
but we can use it in our semiclassical approach to partially restore unitarity. 

The contributions from coefficients f3i with i > [(N + l)/2] can be expressed in terms of coefficients /3, with 
i < [(N + l)/2]. ([x] is the integer part of x.) So, if N is even: 

P(k) = n(k) +detT fj(k), 

N/2-1 Q 

V(k)= £ (3 3 (k) + ^N/2, detT='^. (6) 

If N is odd: 



P(k) =?7(fc)-detT fj(k), 

(N-l)/2 

„(*)= £ m, detT=-f™. (7) 

As a consequence of the imposition of this symmetry on the operator T we obtain two advantages: only traces up 
to half the Heisenberg time tji = N are needed and the eigenvalues are constrained to lie on the unit circle or in 
symmetric pairs. 

These formulae relate P(k) with the traces of powers of T(fc) which, in turn are related semiclassically to periodic 
orbits and to the smoothed density of states |2j|. They have been tested extensively for the hyperbola billiard by 
Keating and Sieber S. 
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B. Green function 



To extend these methods to the calculation of eigenfunctions, we define a generalized Green function G(fc): 

G « = r^w < 8 > 

This operator has poles at the billiard eigenvalues k = fc„ and its residues are the projectors onto the corresponding 
eigenfunctions. It has a Fredholm expression as 

G(k) _ T( t )g(l-T( t )) i (g) 

where C*(l — T(k)) is the transpose of the cofactor matrix of 1 — T(fc) and, as in Eq. (Q), has an expansion in powers 
of T(fc). 

It is then convenient to define a normalized Green operator as 

- Mrr (10) 

where the singularities in the denominator have been eliminated. The normalized Green operator has the property 
g(fc t/ ) = IVvXVvli where \ip v ) is the eigenvector corresponding to eigenvalue k v . Then, we can write g(fc) in the 
following way: 

S (k) = T(fc) C(l-T(fc)) 

SV ; tr(T(fc) C*(l -T(fc))). V ' 

As the cofactor matrix can be expanded in powers of the propagator and as the propagator itself is unitary we write 
the normalized Green operator in terms of the powers of the propagator and their traces up to N/2 (if N is even): 

_V _i V _1 

g(fc) = fiLo 1 c»(fc)T'+ 1 (fc)-dctT(fc) EjLo CiWTt <(fc) 

E£o 1 Ci(fc)tr(T*+i(*)) " detT(fc) eSo' ^(*)tr(Tt *(*)) ' 
where the coefficients Cj(fc) are given by 

f-i 

Ci(fc) = (13) 

An analogous formula can be derived in case N is odd. 

The coefficients Cj(fc) are dependent on the traces of T™ through Eqs. (^) and (13) and thus are independent 



of the chosen representation. On the other hand, the expression for the powers of the propagator will depend on 
the representation chosen for the calculation of the eigenfunctions. If the coordinate representation \q) is chosen, 



Eq. (12) relates the probability density \<j) v (q)\ 2 to the diagonal powers of the propagator (q\T n (k)\q). If the Weyl 
representation is chosen, then Eq. ( |l2] ) gives the Wigner distribution of \<j> v ) in terms of the Weyl propagator [^[Tq. 
Here we choose the coherent state representation to find the equivalent Husimi distribution. 

We remark that Eq. ( |l2[) is a very compact and representation independent derivation of formulae that were 
previously very laboriously derived for the Wigner case. It prepares in an optimal way the grounds for the semiclassical 
approximation because its ingredients are all dependent on classical elements, namely periodic orbits, phase space 
volume and generating function. 

Using the fact that at k — k v the normalized Green function is the projector onto the corresponding eigenstate, we 
can obtain the Husimi distribution as 
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v _i JS _i 

<£jg(Mj£) _ 1 E^o c t (k)(z\T^(k)\z) - dctT(fc) EjLo c,(fc)(z|Tt \k)\z) 
(Z|Z> (Z|Z> ££o l Ci{kMT*+i(k)) - detT(fc) E£o c,(fc)tr(Tt <(fc)) 

This scheme was successfuly applied in simple quantum maps po| . 

III. SEMICLASSICAL APPROXIMATION 

Green's theorem allows us to reduce the Schrodinger equation for the billiard with Dirichlct boundary conditions 
to the following linear homogeneous equation for the normal derivative on the border c/>(s) 

cj)(s) = -l j ds' >(s')K(s, s '; fc), (15) 

where the kernel is 

K(s,s / :k) = ^cos^(s) H^m^-r'is'^), (16) 

with k the wave number, ipis) the angle between the normal at s and the line that connects r(s) with r'(s') (see Fig. 
|l|) and Hi the Hankel function of first type and order one. 
We introduce the wave function /x(s) 



with <f)(p) and n(p) the momentum representations of 4>(s) and /i(s). This transformation makes the kernel symmetric 
and turns Eq. ( |15|) to 

fi(s) = <l>T(s',s;k)tJ,(s')ds'. (18) 

The semiclassical theory of the kernel T(s',s;fc) 0] is based on two fundamental properties: T is semiclassically 
unitary and has an effective dimension N(k) = Lk/ir, where L is the length of the billiard. Moreover, the kernel is 
given by the generating function of the classical Birkhoff map (see Eq.(19)). These properties have been extensively 
tested |26j and will be assumed in what follows. 

Thus, we make the semiclassical approximation by taking T as Bogomolny's operator ( JlSj ) and by evaluating all 
integrals by stationary phase approximation. The T operator for convex billiards in the plane, taking its border as 
Poincare section and using Birkhoff coordinates is 



(-Y 


d 2 l(s',s) 


\2iri J 


dsds' 



1/2 

exp 



(ikl(s',s)-i~v), (19) 



where the bounce map generated by l(s's), the arc length between s and s'; v is the Maslov index. The quantization 
condition is det(l - T(fc)) = 0. 

In the semiclassical theory for the spectral determinant 

P(k) = det(l - T(fc)) (20) 

the approximate unitarity of T can be used efficiently to reduce the number of periodic orbits needed for the com- 
putation of the spectrum. Similar manipulations of the Fredholm formulae allow for the same reduction in the 
semiclassical calculation of single eigenfunctions of the billiard. We write down a formula giving the projector on 
single eigenfucntions as a finite sum involving powers of the map T and its traces. 
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A. Semiclassical traces and determinant 



First we need the traces. It is a well known fact that they adopt the following semiclassical expression 29 

[bnU= idetfj -Mr)\{i/2) ex P ^ kl v ~ *W 2 )) ■ ( 21 ) 

PO,n=n p r 1 ^ P'' 

where the sum goes over all the primitive PO's of the billiard with period n p , which must be a divisor of n, Maslov 
index v p , length l p and monodromy matrix M p . The Maslov index can be interpreted geometrically: ~kv p is the angle 
swept by the unstable manifold of M p along the PO. 
The determinant of T can be obtained as JlJ 

[detT(fc)] sc; = (-l) N cxp(2^AA(fc)), (22) 

where Af(k), the number of states between y k, is 

N(k) = —Ak 2 - —Lk, (23) 

with A the area of billiard and L its length. 

These are the semiclassical ingredients needed for the calculation of the spectrum and of the coefficients Cj (k) . 

B. Semiclassical propagator in coherent state representation 

We hrst obtain the coherent state representation for one iteration of the bounce map: 

(z'\T\z) = J ds ds' (z'\s'} (s'\T\s) (s\z). (24) 

That is to say: 

(z'\T\z) = 



k \ 1/2 / k 



2m 





d 2 i 


j dsds' 


dsds' 



1/2 

exp (ik<Z>(s', s)) (25) 



where 



a) = -z' 2 + l -z 2 -^s> 2 - ^z's' + -V - ^zs + l(s>, s), (26) 
2 2 2a z a 2a z a 

and z' = (q'/cr — iap')/y/2 and z = (q/cr — iap)/V2. The most important contributions come from those points s* 
and s'* that make stationary the phase $: 



9* fa * 1 * *V2 dl „ , 

t-( s > s ) = z + ^r( s > s 

as a a as 



The solution to (|27|) satisfying the reality conditions is 

di 



s =q g^{s ,s )=-p 

dl 

s =q — (s ,s ) =p . (28) 
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This is a classical trajectory from (q,p) to (q',p'). Thus, the matrix element (z'|T|z) will be non zero only if z and 
z' are connected by the classical dynamics. Let us call these points z c and z' c and let us calculate the matrix element 
to next order in their neighbourhoods, (z' c + Sz'\T\z c + Sz). To this effect we expand $(s',s) in Eq. ( |25| ) to second 
order and, after some algebra: 



-iSz z' — iSzz c 



2 Z c z c 



+ 



I- 



+ 



—Sz'Ss' + -Sz 2 + — 
a 2 Za z 



Ss 2 



iy/2 



8z5s 



'2 _ -2 , 2\ 

z c z c ' z c) 

ds' 2 ds'ds ds 2 



(29) 



where 8 s — s — q c and 8s' = s' — q' c . We change to new integration variables 8s and 8s' and obtain: 



(z' c + 8z'\T\z c + Sz)^ UL\ ' (JL} 1 exp(- W /2)x J 





d 2 l 


j d8s d8s' 


dsds' 



1/2 



exp (ik$(8s', 8s)) 



(30) 



We now insert Eq. (p9|) in Eq. (p0|). All terms are constant with respect to integration except the last one in square 
brackets. The resulting integral is the coherent state representation, wit h respect to \8z), of the linearized map, whose 
generating function is quadratic, which we have introduced in Ec. (A5): 



(8z'\T\8z) = — = exp 
Is, 



— (~f c 8z' 2 + 28z'Sz + r c Sz 2 ) 
2s c 

where r c y s c are the matrix elements of the linearized map in complex coordinates. Finally we arrive at: 

-k 



(31) 



x exp 



(z' c + 8z'\T\z c 
1 



Sz) « exp 
1 



2 z 2 



-2 i 2 
Z c + z c 



k ( Sz z c + Szz c + -z c z c + -z c z c 



x exp {ikl — i—i/J (8z'\T\Sz) 



(32) 



This result lets us evaluate the matrix element we were looking for, (z\T n \z)/ (z\z), that will be a sum of contributions 
of periodic points z pp of period n in the semiclassical limit: 



(z|T"|z) 

(z\z) 



pp,n 



z pp + 8z\T n \z pp + Sz) 
(z P p + Sz\z pp + Sz) 



(33) 



To obtain the composition (z pp + <5z|T™|z pp + Sz) we use the expression ( |32] ) and the composition rule of Eq. 
Then 



(z\T n \z) 



^ Xpp exp (iklpp - i-Vpp) x exp 



k 



2s 



-r p p8z 2 + 28zSz + r p 



M 2 ) 



kSzSl 



(34) 



where X pp can be calculated by Eq. ( JA7| ) , v pv — n (because of the Dirichlet boundary conditions) and l pp is the 
length of the PO starting from (q P p,p P p). As we can see, the matrix element behaves as a gaussian in the vecinities 
of the periodic point. This allows us to write the Husimi representation of the n-th power of the propagator as a 
sum of contributions from periodic points of period n. Each term of the sum is a gaussian packet in phase space 
whose parameters are related to the monodromy matrix in complex coordinates. A PO composed by n points will 
give n different contributions to this sum, due to the fact that the monodromy matrices at each point differ. However, 
the invariant properties of these matrices are the same and the usual Gutzwiller- Tabor trace formula p9| can be 
recovered by integration. Of course, a periodic point of period n will contribute also to the rn (r natural) powers of 
the propagator. 
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We should remark at this point that the different semiclassical representations of the propagator in terms od the 
corresponding generating function are only semiclassically equivalent and thus can give different results at finite N. 
This is not true for the calculation fo the spectral determinant, whose semiclassical expression in terms of periodic 
orbits is the same in all representations. It is because of this that the different ways of computing cigcnfunctions are 
not equivalent. For the calculation of |</v(s)| 2 the closed (but not necesarily periodic) orbits are needed [Q. For the 
Wignaer function calculation only periodic points are needed but each contribution is extended in phase space. In the 
present formalism we will obtain the Husimi distributions of eigenfunctions in terms of deformed localized gaussians 
centered in the periodic points, constructed solely in terms of classical information. 



Symmetries 

Our system, the stadium billiard, has two discrete spatial symmetries: R x and R y , the two reflections with respect 
of the coordinate axes. These spatial symmetries in the domain reflect in the border and, thus, in the classical and 
quantum map on it. Their action on the Birkhoff coordinates of phase space (q,p) is 



Rx(q,p) -> (L-q,-p), Ry(q,p) 



R x R y (q,p) 



§+<?, +P 



(35) 



In order to have coherent states on the border with correct symmetries we need to project them using R x and R y , 
the unitary representations of the symmetries R^ja;, y) = \ — x, y) and Ry|x, y) — \x, —y). Then we define 



1 + (j x TL a 



1 + <7yR,y 



(36) 



where a x ,a y = ±1 and and R y move the center of the coherent state according to Ec. (pq). 

In this way, the diagonal matrix elements of the propagator in symmetrized coherent state representation are 



(z\z) 



(z\T 



1 + er x R a 



1 + Uy'R-y 



\z) 



1 



4(z\z) 



((z\T\z) + a x (z\TR x \z) + cr y {z\TRy\z) + ^^(zlTR^R^z)). (37) 



We have already calculated (z|T™|z). We still have to calculate the other three contributions, (z|T"R.,,|,z), 
(z|T n R 1/ |£) and («|T n Ra.Rj,|2;). We can conclude using the results we have already obtained that each of them 
will be a sum of gaussians centered in those points z that the dynamics connects with their symmetric partners, 
R x z, R y z, R x R y z, respectively. These points belong to POs whose periods are 2n which are symmetric under the 
operations R x , R y , R x R y , respectively. The increment RSz with respect to Rz {R = R x , R y , R X R V ) is related to the 
increment Sz with respect to z through: 



RSz = tftSz, 



tn = 



-1 if R = R X 
-1 iiR = R y 
1 if R — R x R y . 



(38) 



Thus we arrive at: 
(dT"RU) 



E 

pp,2n 



X pp exp (ikl p 



7T 
l — Vn 



x exp 



2s z 



(—Tppdz 2 + 2tn5z5z + r P pSz 2 ^ — kSzSz 



(39) 



where the sum goes over the periodic points of period 2n that belong to POs symmetric under R. The quantities s ppi 
r PP , l P p, Vpp and X pp are calculated along the trajectory that connect z to Rz, i.e., half PO. 
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IV. SEMICLASSICAL EIGENFUNCTIONS FOR THE STADIUM 



We use the semiclassical approach we introduced above for the stadium billiard. We choose odd-odd symmetries 
(a x = <jy = —1) and a = 2. We have periodic points with desymmetrized period up to 8 (around 800). We have 
used the symbolic dynamics developed by Biham and Kvale p2| to obtain them. The wave number k is related to 
the maximum period used in the expansion (|l4| ) by P(k) = j^k as 0.4 k. In this way we can obtain semiclassical 
approximations of eigenfunctions of wave number k < 20. 

In Fig. [2] we show the phase space representations of the first six powers of Bogomolny's T operator for k = 20; in Fig. 
|3|, the semiclassical approximations. We see that the exact representations show global maxima in the bouncing ball 
region that cannot be reproduced semiclassically for the lowest powers. However, the overall semiclassical behaviour 
is very close to that of the exact representations. (Because of the symmetries we chose, we have no semiclassical 
approximation to the first power of the operator because the contribution of the only periodic point of period 1 is 
zero.) 

We select two energy ranges: k £ [19.1, 20.0] and k £ [20.5, 21.3]. There are 4 eigenenergies in each of these ranges. 
We show the absolute value of the secular determinant, |P(fc)|, for each of them in Figs. Q and |[ The full line is the 
semiclassical approximation, the dashed line is the secular determinant for Bogomolny's operator. The vertical lines 
are the exact quantum k eigenvalues calculated by the scaling method [^5) . We see a good approximation when we use 
the periodic point expansion. The agreement shows that in this region the spectrum is well represented semiclassically. 
(The discontinuities come from the change in dimension of the operator). 

To keep the method consistent we evaluated the semiclassical Husimi expansion in those values of k that minimize 
the semiclassical secular determinant. We see in Figs. || and [?] the exact eigenfunctions (first column) and their 
corresponding semiclassical Husimi representations (second column) obtained as the real part of Eq. ( |T4| ) . The global 
behaviour is well reproduced; however, the finer details are hard to mimic. The bouncing ball region is problematic: 
in some functions (e.g., k = 21.16) some probability leaks to this region. Probably POs with longer periods that 
approximate the bouncing ball orbits could make a better picture for this region. 

One of the advantages of formula (|lj) for the projector is that it has no singularities between eigenvalues. It is 
possible to study continualy its behaviour as a function of k in order to see its sensitivity to changes in k. Some 
properties of the exact distribution as a function of k are: 

• The distribution is positive at the eigenvalues k n . 

• The distribution has N(k) zeros at eigenvalues k n . 

These two properties follow from the fact that H^ v {z,z) is the modulus of an analytic function. The distributions 
between eigenvalues can become negative. In particular, it can be shown that at the value of k that maximizes P(k), 
the distribution is constant. This properties can be used to control the semiclassical approximations. 

In Figs. H and|| we show the behaviour of the distribution (z\g(k)\z) / (z\z) between the semiclassical eigenvalues 
k = 19.18 and k — 19.38. When k — 19.18 the distribution is positive and has well defined minima that approach 
zero. As we move away from the eigenvalue, the distribution changes smoothly. Initialy it moves away from the plane 
g = in the positive direction, then it comes back and turns negative. During this "evolution" it flattens visibly and 
we can't discern its features. At k = 19.32, aproximately the maximum of the secular determinant, see Fig. 0, the 
distribution is constant. At the semiclassical eigenvalue k = 19.38 the distribution is positive again with well defined 
minima. 

We can see from Figs. || and [7] that the semiclassical approximation is relatively good. It is not trivial to obtain a 
positive defined distribution with N zeroes adding several hundreds of gaussians, each with its phase and deformation. 

V. CONCLUSIONS 

Using Fredholm theory we have given a very compact and representation independent derivation of the projector 
on a single eigenfunction for unitary quantum maps. Expressing the projector in the coherent state basis we wrote a 



8 



semiclassical expression for the Husimi distributions of the billiard's eigenfunctions. Each periodic point contributes 
with a gaussian centered in it whose parameters are calculated only with classical information. We should not 
underestimate the difficulties and complexities inherent to this method. Hundreds of gaussian contributions have to 
conspire to make a positive definite distribution with N that approximate the quantum Husimi distributions. 

The projector ([f2]) can be represented in coordinate space. We obtain (q\ip) {ip\q} , whose semiclassical approximation 
can be directly compared to the probability density in the section. This representation has an additional difficulty, 
since the semiclassical approximation is written as a sum over closed trajectories, periodic or not, in configuration 
space. Those that are not periodic are more in number and more difficult to find. Anyway, we can apply our scheme 
for Bogomolny's T(q',(7) operator and compare the results with the exact quantum calculation. In Fig. [l^ we see 
that the approximation is excelent at this level. 

The maximum period P in the expansions is related to the energy in the way P « 0.4fc. Due to the exponential 
proliferation of orbits in chaotic systems, the method cannot be applied for arbitrarily high energies. The measure of 
this proliferation is the topologic entropy W which relates the number Np of POs of a given period P with the period 
itself, N P = exp(WP) @. For the stadium, W « 0.94. Then, for k w 100 we need POs of periods up to P = 40, 
whose number is Np « exp(0.94 x 40) w 10 17 ! 

The Fredholm method we developed is a first step and shows that the eigenfunctions can be described as expansions 
in terms of the periodic points of the underlying classical system. It eliminates the divergencies associated that the 
schemes based on smoothings in energy have. However, the exponential divergence of periodic orbits poses a serious 
practical problem, as discussed in the previous paragraph. This method can only become practical for large fc„ if 
some way of selecting a few "important" orbits at each value of k can be developed. Some results in this direction 
have been obtained by Vergini and Carlo |^ , ^5| . 



APPENDIX A: COMPLEX PHASE SPACE 



We introduce the following symplectic transformation Z, depending upon parameter a, acting on a point of classical 
phase space (q,p) Bp: 



z \ fq\ = ( l/V2a -«r/V2\ / q 

Pz J \P J \-i/V2a a/V2 J\p 

Imposing reality conditions on the inverse transformation we see that z — ip z . A linear transformation M = 
in (q,p) phase space has a representation M z in (z,p z ) phase space by conjugation with Z 

u/ ■ <i) ; (-4 - 



M z = ZMZ~ L = .- with 
ir s 



(Al) 



a b 
c d 



(A2) 



This {z,p z ) phase space allows a passage to quantum mechanics. This is done in a Hilbert Bargmann space by 
introducing operators z and p z that satisfy the conmutator relations 

[z,p z ] = i%, [z,z] = [p z ,p z ] = 0. (A3) 

Any vector in Hilbert space can be represented in this new space as (z\ip) = J dq(z\q)(q\ip), where the coher- 
ent states are (z\q) = (1 /(nha 2 )) 1 / 4 exp (~(l/h) (z 2 /2 + q 2 /(2a 2 ) - \/2zq/o)). The scalar product is (ipi\ifo) 
f ipi(z) ip2(z) d/i(z) with norm dfi(z) = — exp (— zz/K) dRe(z) dlm(z). 
We now define the Husimi representation of a vector ip as 

_ |(z|^)| 2 
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It is a real positive function for every z in the complex plane. 

The representation of a linear symplectic transformation M in phase space in terms of a unitary operator of Hilbert 
Bargmann space is p3| 



M 



a b 
c d 



(z'\U(M)\z) 



1 



exp — arg(s) I exp 



2s 



{-fz 12 + 2z'z + rz 2 ) 



This representation is up to a phase |2^J2J] and its composition law is 

(z'\\J{M 1 )\z)(z\V{M 2 )\z") dft(z) = X(M 1 ,M 2 ,M 1 M 2 ) (z'\\J(M 1 M 2 )\z"), 



with 



A(Mi, M 2 , MiM 2 ) = exp 



arg(s) - arg(si) - arg(s 2 ) - arg 



SlS 2 



± 1. 



(A5) 



(A6) 



(A7) 



The accumulated phase due to succesive transformations leads to the Maslov index of the trajectory. 

In case the phase space shows periodicity in coordinate or momentum, we have to periodize the coherent states as 

in 
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APPENDIX: FIGURE CAPTIONS 



FIG. 1. Geometry of the billiard. 

FIG 2. Phase space representations for the first six powers of Bogomolny's operator. For each power we show modulus 
(rows 1 and 2) and phase (rows 3 and 4). 



FIG 3. Semiclassical phase space representations for the first six powers of Bogomolny's operator. For each power we show 
modulus (rows 1 and 2) and phase (rows 3 and 4). 



FIG. 4. Secular determinant. In full line we show the PO's approximation; in dashed line, the exact by using Bogomolny's 
operator. For the semiclassical approximation we summed up to period 8. 

FIG. 5. Husimi representations of stadium eigenfunctions (left panels) and their semiclassical approximations (right panels) 
for the energy range of Fig. ^. For the semiclassical approximation we summed up to period 8. 

FIG 6. Secular determinant. In full line we show the PO's approximation; in dashed line, the exact by using Bogomolny's 
operator. For the semiclassical approximation we summed up to period 8. 



FIG 7. Husimi representations of stadium eigenfunctions (left panels) and their semiclassical approximations (right panels) 
for the energy range of Fig. ^. For the semiclassical approximation we summed up to period 8. 

FIG 8. Variation of the normalized Green function in coherent state representation between k = 19.18 and k = 19.28. 
k — 19.18 is a semiclassical eigenvalue. 

FIG. 9. Continuation of Fig. [s| Variation of the normalized Green function in coherent state representation between 
k — 19.28 and k — 19.38. k = 19.38 is a semiclassical eigenvalue. 

FIG 10. Exact eigenfunction (dashed line) for k = 20.5289 and its approximation by using Bogomolny's operator in the 
normalized Green function (full line). 
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